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ABSTRACT 

The correct amplitude and phase modulation formalism of the Blazhko modulation is 
given. The harmonic order dependent amplitude and phase modulation form is equiv- 
alent with the Fourier decomposition of multiplets. The amplitude and phase mod- 
ulation formalism used in electronic transmission technique as introduced by Benko, 
Szabo and Paparo (2011, MNRAS 417, 974) for Blazhko stars oversimplifies the am- 
plitude and phase modulation functions thus it does not describe the light variation 
in full detail. 

The results of the different formalisms are compared and documented by fitting 
the light curve of a real Blazhko star, CM UMa. 
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1 INTRODUCTION 

The periodic modulation of the light curve of a large per- 
centage of RR Lyrae stars, the so-called Blazhko effect, is a 
hundred-year old enigma. The height and time of maximum 
light of these stars oscillate with the same per i od of sev- 
eral days, weeks, months or even years. iBlazhkol (1 19071 ) was 
the first who noticed that no constant period could satisfy 
the observed times of maximum light of RW Dra, an RR 
Lyrae type stai0, and an oscillation in the fundamental pe- 
riod with a cycle length of 41.6 days had to be postulated. 
The striking changes in the height of the light maximum 
and in the shape of the light curve of RR Lyr, whic h turned 
out to be periodic with 40 days, was discovered bv lShaplevI 
(| 19161 ) . Following these discoveries, the Blazhko effect has 
been detected and studied in a number of RR Lyrae stars. 
The large-scale surveys as MACHO and OGLE 



(|Alcock et all 120051 : ISoszvhski et al.l 12001 l201ll) detected 
10-30 per cent incidence rate of the modulation among 
RRab stars in the Magellanic Clouds and the Galactic bulge, 
while a recent groun d-based multicolour photometric survey 
(|jurcsik et al.l I2009T ) and the highly accurat e observations 
of the C0R0T and K epler space missions (jChadidl [20091 : 
iKolenberg et all l201Ch revealed that a significantly larger 
fraction of the fundamental-mode RR Lyrae stars (about 50 
per cent) shows the effect. 

The high occurrence rate of modulated RR Lyrae stars 



1 At the beginning of the 20 th century the origin of the light 
variation of RR Lyrae (cluster-type) variables was unknown and 
they were classified according to the shape of their light curves 
as 'antalgol' type. 



makes the Blazhko effect an even more intriguing problem 
of the pulsation theory. No wonder that it has captivated 
the researchers' interest again in recent years. In spite of 
the fact that the highly accurate space data and the ground- 
based multicolour photometries have given new i nsights into 
the B lazhko phenomenon (for a recent review see lKolenberd 
l201ll ). no generally accepted theory exist s that is able to 
explain the observed features of the effect (|Kovacsll2009l ). 

The connection between the shock-waves in the atmo- 
sphere and light-curve modulation is particularly interest- 
ing. The shock-waves play a significant role in the appear- 
ance of the bump /hump features on the light curve of RR 
Lyrae stars. iPreston. Smak fc Paczviiskil (|l965l ). investigat- 
ing the spectral-feature variations during the 41-day mod- 
ulation cycle of RR Lyrae, suggested a model in which 'a 
critical level of shock-wave for mation moves up a nd do wn' in 
the star's atmosphere. Later, IChadid fc Gilletl |l997l) have 
found that the amplitude of shock-waves is greatly corre- 
lated with the Blazhko modulation. 

To describe the light variations of Blazhko stars, mathe- 
matical models have been propounded. A model that aimed 
to describe th e full Fourier spectrum of th e modulation was 
suggested by iBreger fc Kolenberd (|2006l ) . The triplet was 
interpreted as the non-linear coupling of the main pulsation 
mode and a non-radial mode close to it, and many of their 
combination terms. However, no detailed confrontation of 
the possible predictions of this model with the observations 
(e.g. on the amplitudes of the components of the multiplets) 
has been performed. In other appro aches (|Szeidl fc Jurcsikl 
l2009l : iBenko, Szabo fc Paparol 1201 ll . hereafter BSP11) the 
multiplet structure is the simple result of the modulation of 
purely radial pulsation. 
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In their comprehensive study, BSP11 presented an ana- 
lytical formalism (applied in electronic signal transmission) 
for the description of the light curves of Blazhko RR Lyrae 
stars. Their model shows several light-curve characteristics 
similar to those observed in real Blazhko stars. It was also 
claimed that the new light-curve solution drastically reduced 
the number of necessary parameters compared to the tradi- 
tional methods. However, in BSP11, the suggested method 
was not teste d on real observational data . 

Recently. iGuggenberger et al. I (|2012T ) made an attempt 
at applying the proposed new analytic modulation formal- 
ism to the Kepler data of the complex Blazhko star KIC 
6186029 = V445 Lyr. Although the model light curve showed 
the global properties of the observed one, the fitted light 
curve deviated from the observed data significantly in cer- 
tain phases of the pulsation and the modulation (see figure 
16 in lGuggenberger et al. II2012T ). The surprisingly high vari- 
ance of the residuals was explained by method-specific and 
object-specific reasons. The method did not describe the mi- 
gration of the bumps and humps, and as any 'stationary' 
model, it was unable to follow an irregul ar, time-dependent 
phenomenon (|Guggenberger et al. II2012T ). 

In this paper we look into the possible short-comings of 
the formalism introduced in BSP11. The mathematical for- 
malisms are given in Sect. 2, and the results are documented 
by the different fits of the light curve of a Blazhko star, CM 
UMa in Sect. 3. As the modulation of CM UMa is quite sim- 
ple and regular; no secondary modulation is detected, and 
no higher order than quintuplet modulation components are 
present in the Fourier spectrum, therefore no bias of the 
results arises from any time-dependent irregularity of the 
modulation, what was the case in V445 Lyr. 



2 MATHEMATICAL CONSIDERATIONS 

The Fourier decomposition of the light curve of Blazhko 
RR Lyrae stars is a common technique to analyse their mod- 
ulation properties. If fo and / m indicate the fundamental 
and modulation frequencies, in the typical Fourier represen- 
tation of the light curve, the harmonics (ifo) and the mul- 
tiplets (triplets, quintuplets, septuplets, etc.: ifo ± jfm), as 
well as the harmonics of the Blazhko frequency (fc/ m ) appear 
in the spectrum. 

Let uj = 2nfo and Q — 27r/ m , then the Fourier repre- 
sentation of the modulated light curve is: 



m(t) = m + 2^ b k sin (kQ,t + ip hk ) + 
fc=i 

1+ 

71 r i 

a; sin (iut + ifi) + a J sin (iu)t + jQ,t + ip^j) + 
i=i L 3=1 



<V sin (iu)t - j'Slt + ip ijr ) 



(1) 



where the amplitudes (b, a) and angles (i^b, ip) are constants 
(latter ones are epoch dependent). In the expression of m(t), 
the first sum X]L=i ^fc sin (kftt + ip^k) corresponds to the 
mean light-curve variation during the Blazhko cycle. 



The application of formula 
Eq. [T] leads to the following form 



2] (see Appendix A) to 



m(t) — m + ^2 b k sin (kQt + ipbk) + 
k=l 

n 

J2i a * + / Ai sin l iu)t + <Pi + fvi(t)} ( 2 ) 
1=1 

where a, b, tp and <fib are amplitude and phase constants, 
while f\i{i) and /fi(£) are the amplitude and angle (phase 
or frequency, see Appendix B) modulation functions of 
the ith harmonics of the pulsation. They depend only on 
Q(t) = 27r/ m t and the constant parameters of the Fourier 
series Eq.[T] Therefore, they are periodic functions according 
t° /m 1 (= 27r£l~ 1 ), consequently, they can be approximated 
by the Fourier series: 



fAi(t) = Y2 a % sin (jtit + (p. 



(3) 



and 



/f«(*) = a S' sin + f'H' 

j'=i 



(4) 



where a A , a F , ip A and ip F are amplitude and phase constants. 

In Eqs. [TJ [2] E]and|4j the harmonic orders of the pulsa- 
tion and the different modulation series 71, /, , Z^, are 
limited by the accuracy of the observations. 

The equality of Eqs. [1] and [2] demands that the time- 
history of a modulated star must be equally well described 
by either of them, with about the same variance of the resid- 
ual, no matter which of the two formalisms is chosen. Then, 
the question is raised why does the formalism of BSP11 
give an inferior fit to the Fourier decomposition? If Eq. [5] 
is compared to Eq. 49 of BSP11, the difference is at once 
conspicuous. In Eq. [2] the amplitude and angle modula- 
tion functions (Eqs. [3] and 3]) depend on the harmonic order 
(i) of the Fourier sequence of the unmodulated light curve 
(carrier wave), whereas BSP11 assume that the modulation 
functions depend on the harmonic order in a very restricted 
way. To be more specific, in their approach, each harmonic 
term of the Fourier series of the unmodulated light curve is 
modulated in amplitude and angle by the 



g\{t) = a + 2^ a P > sin f 2-Kp'fmt + (p p > 



p'=i 



(5) 



and 



1 . 
gvi{i) = i ^ a-p sin \ 2Trpf m t + ipp 



(6) 



P =i 



formulae, respectively. We note here that in BSP11, the con- 
stant term of the amplitude modulation function g\(t) is 
identified with the difference between the magnitude and 
intensity mean values, and the constant term of the phase 
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modulation function g?(t) is contracted into the tpi phase 
constant of the pulsation. 

In Eqs.[5]and[6] the BSP11 notation holds, with the ex- 
ception of denoting the harmonic order of the pulsation with 
i, instead of j. Here g\ does not depend on the harmonic 
order of the unmodulated light curve in the least, which is 
inconsistent with Eq. [3] while in gp i , the harmonic order is 
not taken into account correctly with the sole multiplication 
of the modulation function by i (cf. Eq. 

2.1 Consequences 

Due to the adoption of simplified formulae for the amplitude 
and angle modulations, BSP11 drastically reduce the num- 
ber of parameters in their formalism. However, the outcome 
is an inferior fit compared to the Fourier fit (see details in 
Sect. 1 3 - l~f> - The modulation assumed by BSP11 can be effec- 
tuated in electronic signal transmission, but, a priori, it is 
not known how the star modulates its pulsation light curve 
(the 'signal'). 

Another serious problem is the change of the mean 
brightness during the Blazhko cycle. Obviously, the an- 
gle modulation cannot contribute to it. BSP11 state that 
'the formalism naturally explains the mean brightness vari- 
ations', however, in their description it is directly propor- 
tional to the amplitude modulation. Like their amplitude 
modulation formalism describes only an averaged harmonic 
order dependence, it can happen that the mean varies in 
a more complex way than the BSP11 formalism demands. 
In real stars, it is actually observed: e.g., neither the mag- 
nitude nor the intensity averaged variation of the mean V 
brightness of MW Lyr is linearly proportional to the am- 
plitude var iation of the star du ring the Blazhko cycle (see 
figure 14 in lJurcsik et ai]|2008bl ). Moreover, their difference 
is Blazhko phase dependent, thus it cannot be taken into 
account as a sole constant, as in BSP11. 

In principle, the summation of j and j' in Eqs. [3] and 
[4] can start from zero; the j = j' = terms modify the 
amplitude and phase of the unmodulated light curve in each 
harmonic order. These constant terms are included in the 
and parameters of Eq. [2] Therefore, the choice of the 
starting index of the summation from 1 does not have any 
consequence on the results. These extra terms express the 
differences between the mean light curve of a Blazhko star 
and the light curve of an unmodulated RRab star with the 
same physical properties. 

The difference between the phase and amplitude mod- 
ulation functions of BSP11 and the present paper (modula- 
tion functions deduced from the Fourier series) has another 
consequence, too. The observations prove that the amplitude 
versus phase diagrams for the different harmonic orders of 
the light curve during the Blazhko cycle show intrinsically 
different forms (see details in Sect. 13. 2|) . These diagrams 
determined according to the formalism of BSP11 are mor- 
phologically exactly the same for each harmonic order. This 
fact also refers to that the reduction of the degree of freedom 
(decreasing the number of free parameters) by BSP11 alters 
the real phase and amplitude relations of the observations. 

A further important question, which is generally re- 
garded as trivial, is the definition/determination of the 
phase difference between the angle and the amplitude mod- 
ulations. In Eqs. [3] and |4l the angles depend on the starting 
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Figure 1. The V light curve of CM UMa and the light curve 
phased with the pulsation (0.589124 d) and modulation (27.77 d) 
periods are shown in panels A, B and C, respectively. 



epoch, but their difference <&ij = <pfj — ipfj is epoch indepen- 
dent. Since the <&ij values can be different for the possible 
combinations of i and j, the definition of the phase differ- 
ence between the angle and the amplitude modulations is 
not at all obvious. 



3 A REAL CASE STUDY: CM UMA 

In this section, we compare the results of fitting the light 
curve of a real Blazhko star with 

- the Fourier series of Blazhko-frequency-spaced multiplets 
around each pulsation frequency component (Case A, Eq.Q]); 

- the amplitude and phase modulated sum of the pulsation 
light curve (Case B, Eq. [2]); 

- the formalism introduced by BSP11 [Eq. 49 therein] (Case 
C, Eq.EJ. 

In order to index and parametrize the equations in a 
homogeneous manner, Eq. 49 of BSP11 is rewritten in the 
following form: 



m(t) = mo' + 



ao + ^ m sin 

i=i 



A 
a 



af sin (^jflt + tpfj 



2=1 



iujt + i y | ay sin yj'0,t + yJ/J + i 
j'=i 



(7) 



The only difference between Eq. 49 of BSP11 and Eq.Eis 
an additional term: mo'. Without adding this extra constant 
to the equation, a real light curve cannot be fitted. We also 
note that the magnitude zeropoint of Eq. [7] is not identical 
with the zero points of Eqs. [T] and [2] mo' = mo — x ao. 
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CM UMa, one of the Blazhko stars, which was exten- 
sively o bserved in th e course of the Konkoly Blazhko Sur- 
vey II (jSodorl 120121 ) has been chosen for the case study. 
CM UMa was observed with the 60-cm automatic telescope 
of the Konkoly Observatory (Budapest, Svabhegy) in the 
2009 and 2010 observing seasons. The detailed analysis of 
the multicolour light curve of CM UMa will be published 
elsewhere. 

The light curve of CM UMa can be described simply 
with the combination terms of one pulsation (/o) and one 
modulation (/ m ) frequency components. No secondary mod- 
ulation or additional frequency has been detected in the 
Fourier spectrum, which makes CM UMa an ideal target 
for such an investigation. Its light-curve solution consists 
of the Fourier sum of quintuplets around the pulsation fre- 
quency and its harmonics, and the / m and 2/ m modulation 
frequencies. 

The V light curve of CM UMa (shown in Fig. [JJl is 
fitted with different order series of the pulsation (n) and the 
modulation (I, if ,l^l A ,lf ,l A ,l F ) components according to 
Eqs. [JJ [5] and [7] In each solution, the harmonic orders of 
the '+' and ' — ' side multiplets or the amplitude and phase 
modulations are taken to be identical and the same for each 
harmonic of the pulsation, i.e. N = if = l~ = l A = if ; i = 
1, ..rt (Eqs. [1] [3] and 0}. When applying Eq. [71 the orders of 
the amplitude- and phase-modulation series are also taken to 
be identical: N = l A = Z F . The zero point of the magnitudes 
(mo, mo') are also fitted in all the solutions. A low-order 
(I = 1,2) sum of the modulation frequency in Case A and 
B has been taken into account, too. The fits thus comprise 
p = n(4iV + 2) + 21 + 1 parameters when using Eqs. [JJand 
II while p = 2n + AN + 3 for Eq. [JJ 

The pulsation and modulation frequencies are not fitted 
in any of the procedures; fo — 1.697434 and / m = 0.036006, 
determined from the Fourier analysis of the light curve are 
accepted and taken as constants. 

3.1 Comparison of the accuracies of the fits 

Table [JJ compares the rms of the residuals of several light- 
curve solutions using Eqs. [JJ [2] and [7] with different parame- 
ter combinations. In Case A and B, the solutions that have 
the same number of parameters have similar accuracies; the 
rms of the Fourier solution (Eq. [JJ is only marginally, 3 — 5 
per cent smaller than the rms of the amplitude and phase 
modulation formalism according to Eq. [5] This is not at all 
the case, however, when using Eq. [7j Even if similar num- 
ber of parameters are fitted (57/55, 39/43), the residuals 
of Eq. [7] are 60 — 90 per cent larger than the residuals of 
the other solutions. This situation cannot be improved by 
increasing the number of parameters when fitting Eq. [7] as 
the rms of the fit decreases only marginally (by 0.5 per cent) 
when the orders of the pulsation and the modulation are in- 
creased from (16, 2) to (18, 4). If even higher orders (20, 5) 
of the pulsation and the modulation are fitted the decrement 
of the residual is only 0.3 per cent. We note here, however, 
that fitting fifth-order amplitude and phase modulations is 
quite unrealistic, as even in the case of the very complex 
modulation of V445 Lyr, only a third-order FM (frequency 
modulation) a nd a first order AM (ampli tude modulation) 
were assumed ijGuggenberger et al. II2012T ). Also as neither 
pulsation nor modulation components have been detected in 
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Figure 2. The light curves of CM UMa in the lowest and highest 
amplitude phases of the Blazhko modulation (empty circles) and 
the fits according to Eqs. [T] (p=185), H (p=185) , and[7J(p=55) 
are denoted by +, x and • symbols, respectively (top panels). 
The inserts magnify the light curves and the fits at minimum and 
maximum light. The bottom panels show the residuals of the fits. 



the Fourier spectrum at harmonic orders higher than 18, the 
n=20 solution has no relevance compared to the n=18 one. 
As no secondary modulation of CM UMa has been detected, 
therefore the parameter number when applying the formal- 
ism of Eq. 7 cannot be increased this way either. Therefore, 
the p—55 (n=18, iV=4) solution of the BSP11 formalism is 
accepted and used for the comparison. 

The similar accuracy of the solutions of Eqs. [JJ and [2] 
confirms that these forms are, in fact, equivalent in prac- 
tice. Of course, minor difference between these solutions may 
arise from the different role and significance of the param- 
eters in the two formulae. Therefore, the solutions do not 
necessarily have to have exactly the same accuracies when 
the same number of parameters (n, N, I) axe fitted. 

In the top panels of Fig. [2] the light curve in the lowest 
and highest amplitude phases of the modulation are shown 
and their fits using Eq. [JJ (p=185), Eq. [2] (p=185) and Eq.[7] 
(p=55) are overplotted. (The p=63 solution of Eq.[7]does not 
give any noticeable difference if compared to the results of 
the p=55 one.) While the first two fits are indistinguishable 
from each other and they follow even the smallest features of 
the observed light variations accurately, this is not true for 
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Table 1. Comparison of the rms of the residuals when fitting the light curve of CM UMa according to Eqs. [T] [2] and [7] 
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Figure 3. Residual spectra of the light curve prewhitened for 
the solutions of Eqs. |T] (p=185), (p=185), and [7J (p=55) arc 
shown. Note that the y-scales are the same in the three panels. 
The vertical grid indicates the positions of the pulsation frequency 
and its harmonics. 



the fit obtained using Eq.0 Large, 0.02 — 0.05 mag system- 
atic differences, especially around light minima and maxima, 
appear between this fit and the observations. The bottom 
panels show the differences (fit - observation) for the three 
light-curve solutions. While the first two show nothing but 
noise around zero, in the third case, the fit does not follow 
the light variation correctly in any phase of the pulsation. 

The spectra of the three residuals shown in Fig. [3] also 
document that the differences between the results are sub- 
stantial. While the mean levels of the residual spectra are 
around 0.0003 mag in the first two cases, it is as high as 
0.0009 mag in the third one. Moreover, when using Eq. [7] 
0.006 — 0.009 mag amplitude signals at around the positions 
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Figure 4. The observed and fitted amplitudes vs phases of the 
different harmonic orders of the pulsation light curve during the 
Blazhko cycle are shown. The phase values are divided and the 
amplitudes arc multiplied by the number of the harmonic order 
(j) in each plot for a better visualisation and comparison. Note 
that the same, 0.7-rad phase ranges are shown in each panel, while 
the amplitude ranges are different in the different rows of panels. 
The symbols denoting the observations and the three different 
synthetic fits are the same as in Fig. [2] 



of the pulsation frequency and its harmonics remain in the 
residual. 
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3.2 Comparison of the amplitude and phase 
relations of the fits 

The maximum phase - maximum brightness diagrams of 
Blazhko RR Lyrae stars show a large variety. The di- 
rection of going around these loops and their shape are 
defined by the relative amplitudes and phases of the 
maximum-brightness and maximum-phase variations dur- 
ing the Blazhko cycle. The same is true for the shapes of 
the maximum phase - maximum brightness loops in the dif- 
ferent harmonic orders of the pulsation. These curves can be 
determined from the amplitude and phase variations of the 
different harmonics of the pulsation light curve during the 
Blazhko cycle. 

The amplitude and phase variations have been shown to 
have different shapes in the different harmonic orders even 
in a very regular B lazhko star, MW Lyr (see figure 12 in 
lJurcsik et a"ill2008al ). Striking differences between the vari- 
ations of the amplitudes of /o and 2/o (Ai and R21) of RR 
Lyr h ave been also shown (see figure 12 in iKolenberg et all 
Ell). 

In this section, the amplitude and phase variations in 
the different harmonic orders of the light curve of CM UMa 
are shown, and compared to the same plots derived from 
synthetic data according to the three light-curve solutions 
introduced in the previous section. 

Using the parameters of the p=185, p=185 and p=55 
solutions of the three different formalisms, synthetic light 
curves of CM UMa have been generated. Dividing the syn- 
thetic data into small, homogeneous Blazhko phase bins, the 
variations of the pulsation light curves during the Blazhko 
cycle predicted by the different fits can be followed and com- 
pared to the observations. Both the observed and the fitted 
pulsation light curves in the different Blazhko phases are 
then fitted with a 15th-order Fourier sum of the pulsation 
frequency. The phases and amplitudes of the /o,2/o,..i/o 
components characterize the modulations of the pulsation 
components in the different orders; they define the maxi- 
mum phase-maximum brightness loops. 

In Fig. 31 the amplitude vs phase variations of the first 
nine harmonics of the pulsation during the Blazhko cycle 
are plotted. The plotted data are derived from 12 and 30 
Blazhko-phase bins of the observations and the fits, respec- 
tively. The phase variations are normalized by the division 
by the harmonic order i (cf. Eq. [7} and the amplitude varia- 
tions are multiplied by the harmonic order for better visual- 
isation. The synthetic data using Eqs.[T]and [2]fit the obser- 
vations in the first six orders absolute correctly. In the 7-9th 
orders, minor differences both between these fits and the fits 
and the observed values appear, because the observed values 
are somewhat uncertain in these harmonic orders. 

The observed and fitted (via Eqs. [T] and [2} loops are 
significantly different in the different orders; the most strik- 
ing differences are in the phase ranges of the loops. The 
phase ranges, even if normalized by the division of the or- 
der number vary between 0.1 and 0.6 rad. Large diversity in 
the morphology of the loops are also evident, e.g. it is 'egg'- 
shaped in the first order while in the 2nd and 6th orders the 
loops are degenerated, they are reduced to one-dimensional 
curves in these harmonic orders as the amplitude and phase 
variations (see Fig. [4| are quite symmetrical in these orders, 
and their phase differences are 180° . 
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Figure 5. Variations of the Fourier amplitudes and phases of the 
harmonic orders of the pulsation light curves during the Blazhko 
cycle. From the left to the right, the panels show the results ob- 
tained directly from the observations and from the three different 
fits. Filled symbols denote the phase variations. 



The synthetic data corresponding to Eq.[7]do not follow 
the real variations of the phases and amplitudes of the max- 
ima with the harmonic orders, due to the strong restrictions 
of the variations in the amplitude and phase modulations 
of the formalism of BSP11. The morphology of these loops 
are identical in each order, as differences only in the relative 
scaling of the amplitudes are allowed. 

Fig. [5] documents the phase and amplitude variations 
of the pulsation light curve versus Blazhko phase in the dif- 
ferent harmonic orders. The observations and the results of 
the three fits are plotted in separate columns here. Differ- 
ent amplitude and phase ranges are shown for the different 
harmonic orders but the same ranges are used for the ob- 
served and fitted values in each order (in each row in Fig. [5}. 
Again, we can see that the first two fits give back the ob- 
served variations similarly well, with high precision. In Case 
C, however, neither the mean values nor the amplitudes and 
shapes of the predicted variations match the observations 
correctly. The BSP11 formalism is not capable of fitting the 
observed differences of the phase relations between the am- 
plitude and phase variations in the different orders. 



4 CONCLUSIONS 

The Fourier sequence describing the light curve of a Blazhko 
star (Eq. [1]) has been transformed into the form of an am- 
plitude and angle modulated signal (Eq. The two de- 
scriptions are fully equivalent to each other. The correctly 
deduced amplitude and angle modulation functions, fAi(t) 
and /f;(£) depend upon the harmonics of the unmodulated 
light curve (i). Since these are periodic functions, they can 
be expressed as Fourier sums (Eqs.[3]and|4|). A priori, there is 
no knowledge about the parameters of these equations, they 
can be determined only through observations. The possible 
correlations among them can be revealed by observations, 
as well. 

In their study, BSP11 employed an amplitude and angle 
modulation pattern used in electronic signal transmission. 
(In this technique, the coding of the modulation of a signal is 
known in advance and can be controlled.) They assume that 
the amplitude modulation function does not depend upon 
the harmonic order of the unmodulated signal (light curve) 
and the difference between the angle modulation functions 
in the different orders is simply the multiplication by the 
harmonic-order number. This oversimplified procedure leads 
to the drastically reduced number of parameters but, on the 
cost of an unacceptably poor fit. 

The analysis of the light curve of the Blazhko RR Lyrae 
star, CM UMa (Sect.HJ), clearly shows that the BSP11 for- 
malism gives a poor fit, especially around the minimum and 
maximum light and the ascending branch of the light curve. 
Even if the number of the parameters are increased, the 
fit does not improve. It should be emphasized that these 
parts of the light curve reflect the changing strength and 
occurrence of the shock-waves during the Blazhko cycl e 
l|Preston. Smak fc Paczvhskil Il965l ; IChadid fc Gilletl Il997t ). 
We thus conclude that the BSP11 approach does not take 
the nonlinear interactions, which determine the final form of 
the light variation, fully into account because of the strong 
restrictions of the formalism. Its capability to describe the 
changes of the pulsation light-curve's shape is seriously lim- 
ited, since it only shifts the pulsation light curve in phase 
and scales it in amplitude periodically during the modula- 
tion. 

Plots that show the amplitude and phase variations 
of the different harmonic orders during the Blazhko cycle 
(Figs.|4]and[5]) reveal it convincingly that the amplitude and 
angle modulation functions given in BSP11 are unsuited to 
describe the Blazhko modulation correctly. 

Nevertheless, the Blazhko modulation can be cor- 
rectly interpreted as amplitude- and angle-modulated sig- 
nal (Eqs. [2HU- This fact hints at the possibility that dur- 
ing the Blazhko cycle the radial pulsation of an RR Lyrae 
star is modulated in phase/frequency and amplitude, and 
the fundamental period is subject to real oscillation. Up to 
now, the only model that is in conformation with this sce- 
nario has been suggested by IStothers! (|2006l ), although the 
inside physics of thi s mod el has been strongly critic i zed by 
ISmolec et ail |201ll) and iMolnar. Kollath fe Szabl (|2012T ) 
recently. 
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APPENDIX A: TRIGONOMETRIC 
IDENTITIES 

From the well-known identities of trigonometry, we obtain: 

Ai sin (lj + xi) + A 2 sin (lj + X2) = A sin (w + \) (Al) 
where 

A = \]a\+A\ + 2Ai A 2 cos (X2 - Xi ) 

and 

Ai sinxi + A 2 sinx2 
tan X = 1 : — 

Ai cos xi + A 2 cos \ 2 

By reiterating the formula I All we come to the relation 

71 

^2 At sin (u + X i) = A sin (u + X ) (A2) 
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where both A and \ are bounded functions of 
{Ai, A2, A n ,xi,X2, Xn)- This relation holds even if Ai 
and Xi are time dependent. 



APPENDIX B: FREQUENCY AND PHASE 
MODULATIONS 

Two types of angle modulation exist, frequency (FM) and 
phase (PM) modulations. 

In the case of FM, the frequency of the unmodulated 
wave (carrier signal), fo, is modulated by the modulating 
signal, Qp(t). The instantaneous frequency is 

F(t) = f + k F Q F (t), (Bl) 

where kp constant depends upon the modulating system. 
The instantaneous phase of the modulated wave is 

9(t) = 2nf t + 2nk F [ 6 F (r)dr. (B2) 
Jo 

(For simplicity, zero phase is assumed at t = 0.) 

In the case of PM, the phase of the modulated wave is 

#(t) = 2-Kf t + fc P e P (t), (B3) 

where <dp(t) and kp are the modulating signal and a con- 
stant, respectively. The instantaneous frequency of a PM 
wave is 

F(t) = f + k P ^±. (B4) 

In both cases, the wave's frequency and phase vary from 
moment to moment. From Eqs. IB 1 1 and IB4I follows that the 
two modulations describe the same modulated signal, if the 
following relation holds: 

e,m-£^2. (») 

The only difference between the two descriptions is that in 
the case of FM, the frequency modulation of the carrier is 
given by the time derivative of the PM modulated signal. 
Unless some information is available about the modulation 
in advance, (as it is the case in electronic signal transmis- 
sion), it may not be obvious, which one of the two types (FM 
or PM) is realized. If no information is available regarding 
the angle modulation (as in the case of Blazhko modula- 
tion), it is impossible to identify it as an FM or PM signal. 
This explains why the expression 'angle modulation' is used 
thorough this paper. 



